Uniformly frustrated XY model without vortex-pattern ordering 
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The uniformly frustrated XY model with / = 1/3 on a dice lattice is shown to possess a so well 
developed accidental degeneracy of its ground states that the difference between the free energies 
of fluctuations does not lead to the stabilization of a particular vortex pattern down to zero tem- 
perature. Nonetheless, at low temperatures the system is characterized by a finite helicity modulus 
whose vanishing (at a finite temperature) is related with the dissociation of half-vortex pairs. 

PACS numbers: 74.81. Fa, 64.60.Cn, 05.20.-y 



It is well known that in the presence of external mag- 
netic field a regular two-dimensional Josephson junction 
array can be described W the Hamiltonian of a uniformly 
frustrated XY model [ij: 



H = -JY^ cos(^k - - ^jk) , 

(jk) 



(1) 



where J is the Josephson coupling constant, the fluctuat- 
ing variables ipj are the phases of the order parameter on 
superconducting grains forming the array, the quenched 
variables Ajj. are determined by the vector potential of 
the field and the summation is taken over all pairs of 
grains connected by a junction. The form of Eq. as- 
sumes that the currents in the array are sufficiently small, 
so their proper magnetic fields can be neglected. 

The directed sum of the variables Ajk = — Akj along 
the perimeter of each plaquette should be equal to /, the 
ratio of the magnetic fiux per plaquette to the flux quan- 
tum $0- It is sufficient to consider the interval / S [0, 
because all other values of / can be reduced to this in- 
terval by a simple replacement of variables . The case 
of / = corresponds to the absence of frustration. 

For rational / the nature of the vortex pattern in 
the low-temperature phase of a uniformly frustrated XY 
model is usually unambiguously determined by the struc- 
ture of its ground states. The best known examples be- 
longing to this class are the models with f — 1/2 and 
square or triangular lattice 0. On the other hand, quite 
often the ground states of a uniformly frustrated XY 
model are characterized by an accidental degeneracy not 
related to symmetry In such cases the nature of vor- 
tex ordering at low temperatures cannot be determined 
without comparing the free energies of fluctuations in the 
vicinities of different degenerate ground states 0. 

Quite often (for example, when the lattice is triangular 
and 1/3 or 1/4 .4;]), this mechanism of the removal of 
an accidental degeneracy works already at the harmonic 
level, but in some situations (at / = 1/2 on honeycomb, 
dice and kagome lattices) one has to take into account 
the anharmonicites (6|. However, in all the cases inves- 
tigated insofar, at low enough temperatures a particular 
periodic vortex pattern can be expected to be stabilized 
by fluctuations (in the thermodynamic limit). 



In the present work we demonstrate that the uniformly 
frustrated XY model with / = 1/3 and dice lattice (see 
Fig. 1) has rather unique properties. Namely, this is the 
first example of a frustrated XY model in which the ac- 
cidental degeneracy of ground states is so well developed 
that vortex pattern remains disordered at arbitrarily low 
temperature, although the free energies of fluctuations 
are different for different periodic patterns. As a conse- 
quence, the only phase transition which takes place in 
this system with the increase of temperature is related 
to dissociation of pairs of fractional vortices with topo- 
logical charges ±1/2. A possibility for the stabilization 
of a specific periodic pattern can appear only if one goes 
beyond the limits of the XY model. In the conclusion we 
briefly discuss the removal of the accidental degeneracy 
by magnetic interactions and its possible consequences. 

The interest to magnetically frustrated superconduct- 
ing systems with a dice lattice have been motivated by 
the unusual properties of a single electron spectrum in 
this geometry [3. In recent years superconducting wire 
networks and Josephson junction arrays with a dice lat- 
tice have become the subject of active experimental inves- 
tigations H, 0, 0| and numerical simulations |0| . Mag- 
netically frustrated Josephson junction arrays formed by 
rhombic plaquettes have been also discussed in the con- 
text of creation of topologically protected qubits [l^ . 

Since both ipj and Ajk depend on a particular choice 
of the gauge, it is more convenient to describe different 
states of the system in terms of the gauge-invariant phase 
differences 
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FIG. 1: Dice lattice is periodic and consists of identical rhom- 
bic plaquettes with three different orientations. 
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which below are always assumed to be reduced to the 
interval (— 7r,7r). The subscripts j and k will be used to 
denote six- and threefold coordinated sites respectively. If 
one defines the variable mjj' to be given by the directed 
sum of variables 6 over the perimeter of the plaquette 
(jkj'k') divided by 27r, the plaquettes for which mjj' is 
equal to 1 — / (rather then to — /) are usually referred to 
as containing vortices. Different local minima of can 
be then classified by specifying the positions of vortices, 
whose concentration for / = 1/3 should be exactly equal 
to one third. 

The remarkable property of the considered XY model 
is that all minima of ^ corresponding to different vor- 
tex configurations in which vortices do not occupy ad- 
jacent plaquettes have the same energy (per bond), 
i?o — —(2/3) J. Some of these states are schematically 
shown in Fig. 2. In all such states the variables 9 are 
equal to 7r/3 on all bonds surrounding a vortex and to 
zero on all other bonds (shared by two plaquettes with- 
out vortices). It is evident that these states correspond 
to the absolute minimum of energy, because they mini- 
mize the energy separately for each plaquette containing 
a vortex and for each of the remaining bonds. The bar- 
riers between the "neighboring" ground states have the 
height Ex, = 6(2 - V3) J ~ 1.61J. 

All these ground states [whose degeneracy survives 
even if the interaction in Q| deviates from cosine] can 
be put into correspondence with the ground states of the 
antiferromagnetic Ising model defined on the triangular 
lattice T formed by the sixfold coordinated sites {j}. In 
all ground states of such a model each triangular plaque- 
tte should contain exactly one bond with parallel spins 
(sjSj' = 1) The existence of a mapping becomes 

evident as soon as one notices that in all ground states 
of the considered XY model each hexagon formed by 
three neighboring plaquettes should contain exactly one 




FIG. 2; Some ground states of the frustrated XY model with 
a dice lattice and / = 1/3. The black plaquettes correspond 
to m = 2/3 and white plaquettes to m = —1/3. 



vortex. The plaquettes with vortices can be then identi- 
fied with the bonds connecting parallel spins by setting 
SjSj' = 2mjj, - 1/3. 

At zero temperature, T = 0, the antiferromagnetic 
Ising model on a triangular lattice is characterized by an 
algebraic decay of correlation functions and a finite 
extensive entropy, which is related to the possibility of 
creation of zero-energy domain walls forming closed loops 
[l^ . The same partition function can be interpreted as 
the partition function of the SOS (solid-on-solid) model 
suitable for the description of the height fluctuations on 
the (111) facet of a crystal with a simple cubic lattice [l5| . 
In terms of this SOS model zero-energy domain walls cor- 
respond to zero-energy steps, and the algebraic correla- 
tions of spins are translated into logarithmic correlations 
of integer variables /ij and nk, which can be associated 
with the height of the surface. These variables can be 
introduced following the relations 

ft-j±e„ = ± 3TOjj±e„, = ^ ^ /ij , 

j=j(k) 

where (with a — 1,2,3) are the three basic vectors 
(ei + 62 + 63 = 0) of T and j(k) are the three nearest 
neighbors of k on the dice lattice. According to Ref. 0, 
for |ki - kzl » 1 

((riki -rikj') « -^ln|ki -kal . (2) 

The form of Eq. (0) demonstrates that the SOS model 
is in the rough phase and that at T = the large-scales 
fluctuations of n can be described by a continuous Gaus- 
sian Hamiltonian, 

H = 11 Sv[Vnf , (3) 

where the dimensionless effective rigidity K (which is of 
entropic origin) is equal to Kq = tt/9 The phase 
transition of the SOS model to the smooth phase would 
take place when K ~ 11/2 thus the system is situated 
relatively far from the transition point. 

The simplest periodic ground state shown in Fig. 2(a) 
(the striped state) has been discussed in Ref . IllI In terms 
of the SOS representation this state has the maximal 
possible slope, whereas a flat state of the SOS model 
(with rik = const) corresponds to the honeycomb vortex 
pattern of Fig. 2(b). Fig. 2(d) shows two parallel steps (of 
opposite signs) which separate flat states with An = ±1. 
If the left step of Fig. 2(d) is repeated as often as possible, 
one obtains the striped state of Fig. 2(a). On the other 
hand, the repetition of the right step of Fig. 2(d) leads to 
the zigzag state shown in Fig. 2(c). However, the steps 
do not have to be straight, and a typical ground state 
looks rather disordered, Fig. 2(e). 

If the SOS model would be in a smooth phase, the 
vortex configuration would be characterized by the long- 
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range order corresponding to the formation of the hon- 
eycomb pattern of Fig. 2(b). The fluctuations of rik and 
/ij lead to the replacement of the true long-range order 
by an algebraic decay of correlations of rrijy, 

Hiji™j2j^> « |ji - j2r'' , 

where rj = 2, as follows from the results of Ref. ^3- How- 
ever, one can expect this correlation function to be mod- 
ulated according to the honeycomb pattern. 

At finite temperature, T > 0, the equivalence between 
the considered XY model and the Ising model is no 
longer exact, because the XY model (i) allows for the 
existence of continuous fluctuations (spin waves) and (ii) 
has a more complex classification of topological defects. 
Numerical calculation of the integrals over the Brillouin 
zone which determine the free energies of harmonic fluc- 
tuations in the vicinity of the three periodic ground states 
shown in Figs. 2(a), 2(b) and 2(c) reveals that this free 
energy is the lowest in the honeycomb state The 
difference between the free energies of fluctuations in the 
zigzag and honeycomb states (normalized per single site 
of T) is given by 7T, where 7 « 2.27 • 10"^. 

Since in terms of the SOS model the zigzag state can 
be considered as the sequence of steps with the unit 
density, the same quantity can be also used as an es- 
timate for the effective energy of a step per elementary 
segment, Est ~ jT. The positiveness of Est should lead 
to a decrease of fluctuations of rik, but since one al- 
ways has Est/T <C 1, this decrease has to be relatively 
small. In terms of Q the influence of a small step energy. 
Est ~ jT, is translated into a very small (< 1%) correc- 
tion to Ko, K = Kq + (2/-\/3)7. This definitely leaves the 
system far from the transition to the smooth phase, and, 
therefore, cannot lead to any qualitative changes from 
the zero-temperature behavior. 

In terms of the original XY model each step of the 
SOS model corresponds to a line whose crossing shifts the 
phase variables ipj by tt (with respect to what they would 
be in the absence of this step) on one of the three sublat- 
tice into which the triangular lattice T can be split. After 
crossing three steps of the same sign the variables tpj are 
shifted by tt for all three sublattices. In other words, the 
state which is obtained in such a way differs from the 
original state by a global rotation of all phases by tt. As 




FIG. 3: A possible structure of a half-vortex. 



a consequence, the elementary topological excitation of 
the considered XY model is an object where three steps 
of the same sign merge together (see Fig. 3) that looks 
like a screw dislocation with Burgers number b = ±3, on 
going around which the phase experiences a continuous 
rotation by tt. For brevity we shall call such defects half- 
vortices. Since the phase shift by tt can be achieved by 
a phase rotation in both directions, the signs of the two 
topological charges of a half-vortex (the Burgers number 
and the vorticity) are not related to each other and can 
be arbitrary. 

The core of a half-vortex can be associated with the 
three-plaquette cluster which instead of containing ex- 
actly one vortex contains either two vortices or no vor- 
tices at all (see Fig. 3). In the framework of the antiferro- 
magnetic Ising model the analogous defect cannot exist, 
because each plaquette can contain only an odd number 
(1 or 3) of bonds with s^sy = 1. The core energy of a 
half- vortex should be of the order of J. 

In accordance with the double nature of half-vortices 
their interaction consists of two contributions of differ- 
ent origin. Both of them are logarithmic. The first one 
(the direct interaction) is related to the energy which 
is required to create the phase twist around the cores 
and is completely analogous to the interaction of ordi- 
nary vortices (with integer vorticity) . This interaction is 
characterized by the prelogarithmic factor Py = (7r/2)r, 
where F is the helicity modulus of the system [at T = 
in the honeycomb state F = Fg = (5/4-\/3)J]. The sec- 
ond contribution is of entropic origin and is related to 
the interaction of half-vortices as dislocations. It follows 
from Eq. (PJ that for this interaction the prelogarithmic 
factor is given by 

Kb'^ 

Puib) = —T, (4) 

where one should put 6 = 3. 

At low temperatures (T <C J) the direct interaction 
of half- vortices is dominant, which binds them into small 
pairs with zero total vorticity. However, the Burgers num- 
ber of such a pair does not have to be zero, but can be 
also equal to ±6. This returns one to the situation in the 
antiferromagnetic Ising model in which the elementary 
topolog ical excitations are the dislocations with b = ±6 
|l6l Il8j , the only difference being a slightly larger value 
of K. Substitution of X « iCo = 7r/9 and 6 = 6 in Eq. Q 
gives fb(6) ~ 2, which is insufficient for such dislocations 
to be bound in pairs The application of the Debye- 
Hiickel approximation to the two-dimensional Coulomb 
gas formed by dislocations shows that cd , the concentra- 
tion of free dislocations, should be exponentially small in 
l/T. The value of cd determines a temperature depen- 

— 1/2 

dent correlation radius fc(T) oc , beyond which K is 
renormalized to zero and even the algebraic correlations 
of vortex pattern are destroyed. 
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The presence of free dislocations leads to the screen- 
ing of the entropic part of the logartihmic interaction of 
half-vortices at the scales which are large in comparison 
with rc{T). Nonetheless, at low temperatures the sys- 
tem will be characterized by a finite value of r(T), since 
all half-vortices will be bound in pairs by their direct 
interaction. With increasing temperature a phase transi- 
tion will occur related to appearance of free half- vortices 
and vanishing of r(T). It will be completely analogous 
to the Berezinskii-Kosterlitz-Thouless phase transition in 
the conventional XY model (without frustration), the 
main difference being that half-vortex pairs dissociate 
when T = (7r/8)r(r) jl^, whereas the pairs of ordinary 
vortices dissociate when T = (7r/2)r(T) 0. 

Thus, in the present work we have shown that in the 
frustrated XY model with a dice lattice and / = 1/3 the 
vortex pattern is disordered at any temperature (becom- 
ing quasi-ordered only at T — 0). Nonetheless, at low 
temperatures the helicity modulus is finite and jumps to 
zero only at T = Thv ~ (7r/8)ro « 0.28 J, where the 
pairs of half-vortices dissociate. This estimate is not far 
from the value of the transition temperature, Tc ~ 0.2 J, 
obtained in numerical simulations of Ref. ilL 

A possibility of vortex pattern ordering appears only 
when one goes beyond the limits of the XY model and 
takes into account some additional mechanism of the re- 
moval of an accidental degeneracy. In the case of a prox- 
imity coupled array the main role will belong to the 
energy related to the magnetic fields of currents in the 
array, which is minimized in the striped state |2ll |. In 
terms of the SOS model this leads to a negative step en- 
ergy. Est < 0. In the limit of weak screening, when the 
corrections to currents from their proper magnetic fields 
can be neglected, E^t = —j^J'^/E^, where E'* = ^l/A-K^a 
is the characteristic energy, a is the lattice constant (of 
a dice lattice) and the numerical coefficient fi for both 
types of steps shown in Fig. 2(d) can be written as 



H « sin^(7r/3) • [A2 - A4 - 2(A5 - Ag 



0.25 



Xi = —Li/a > being the dimensionless values of mutual 
inductances, L,;, between dice lattice plaquettes jl^l- For 
a KiSum ^ E^!v lO'^K. 

With decrease of T the ratio \ Est\/T is increased, which 
for \Est\/T <C 1 will manifest itself only in the decrease 
of K in Eq. With further decrease of T a phase 
transition can be expected to occur to a phase with a 
non-zero slope. In terms of vortices this phase will be 
characterized by a true long-range order manifesting it- 
self in deviation of occupation probabilities for plaquettes 
with different orientations from 1/3. This phase transi- 
tion can be expected to happen when |-Est| T, i.e., 
T/J - ifJ.T/E^)^/'^ - 10^2^ and, according to the re- 
sults of Ref. I23, has to be of the first order. However, 
at T/J ^ 0.05 the relaxation of vortex pattern is likely 
to be dynamically quenched [at T/J — 0.05 one has 



exp(— i?b/T) 10 ], which may prevent the observa- 
tion of such an ordering in experiments or simulations. 
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